Optimal Time-Reversal of Multi-phase Equatorial States 
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Even though the time-reversal is unphysical (it corresponds to the complex conjugation of the 
density matrix), for some restricted set of states it can be achieved unitarily, typically when there 
is a common de-phasing in a n-level system. However, in the presence of multiple phases (i. e. a 
different de-phasing for each element of an orthogonal basis occurs) the time reversal is no longer 
physically possible. In this paper we derive the channel which optimally approaches in fidelity 
the time-reversal of multi-phase equatorial states in arbitrary (finite) dimension. We show that, 
in contrast to the customary case of the Universal-NOT on qubits (or the universal conjugation in 
arbitrary dimension), the optimal phase covariant time-reversal for equatorial states is a nonclassical 
channel, which cannot be achieved via a measurement/preparation procedure. Unitary realizations 
of the optimal time-reversal channel are given with minimal ancillary dimension, exploiting the 
simplex structure of the optimal maps. 
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I. INTRODUCTION 

> . 

V^Q , Time reversal is not a physically achievable transformation on arbitrary quantum states, since it corresponds to a 
positive, but not completely positive map [lj. However, if restricted to special sets of states, time reversal can be easily 
achieved by a unitary transformation, e. g. when generating NMR spin-echoes by turning over the fixed magnetic 
field along the rotation axis, with all the spins rotating in the equatorial plane. For a general n-level system, this is 
possible only when the levels are equally spaced, corresponding to a common de-phasing. However, more generally, 
when the levels are not equally spaced, we have different de-phasing for each level, and in the presence of multiple 
phases the time reversal is no longer physically possible. 

The ideal time-reversal of a generic state corresponds to the complex-conjugation (or, equivalently, transposition) 
of the corresponding density matrix. Such transformation has also recently attracted much interest in relation to 
the problem of entanglement, in regards to the so-called PPT (positive partial transpose) criterion 0, 0- Since 
complex-conjugation cannot be achieved unitarily, one can try to approximate the transformation with a physical 
channel, optimizing the fidelity of the output state with the complex-conjugated input. For the set of all pure states 
Q" 1 the resulting optimal channel is " classical" 0, [B| , in the sense that it can be achieved by state-estimation followed 
by state-preparation. In this paper we show that for multi-phase equatorial states the optimal phase covariant time- 
• i-h , reversal for equatorial states is a nonclassical channel, namely it cannot be achieved via the measurement/preparation 
procedure. We will see that the optimal channels form a simplex (i. e. a convex set which is generated by convex 
combination of a finite set of extremal points, e. g. a tetrahedron). Such a structure simplifies the search for unitary 
realizations of the channels, which will be derived in the following for minimal ancillary dimension. 

The paper is organized as follows. In Section [n] we introduce the notation, and derive the optimal multi-phase 
conjugation maps and the corresponding fidelity. In Section lTTTl we compare the present optimal phase-covariant maps 
with the universally covariant ones, and discuss their relation with optimal state-estimation and phase-estimation. 
In Section IIVI we analyze the simplex structure of the set of optimal multi-phase conjugation maps, and explicitly 
construct their unitary realizations with minimal ancilla dimension. Section [3 closes the paper with some concluding 
remarks. 
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II. OPTIMAL MULTI-PHASE CONJUGATION MAPS 



In the following we will restrict attention to equatorial states of a d-dimensional quantum system, defined as 

m^})) = i(|0) + e** |1) + e* |2> + ... + e^- |d - 1)), (1) 

expanded with respect to the fixed orthonormal basis 

B = {|0),|l),...,|d-1)}. (2) 

of the Hilbert space H of the quantum system. We consider transformations that treat all input states Q in the 
same way, namely that are covariant under the group U(l) x ' d_1 ) of rotations of d — 1 independent phases {4>j}. The 
state Q can be equivalently written as 

W{&})> = UdhjMo), (3) 

where 

d— 1 

\i>o)=d- 1 / 2 J2\i) (4) 

i=0 

is a fixed real state and 

U({4> j }) = \0)(0\+Y i e^m\ (5) 
i=i 

is the generic element of U(l) x ^ _1 - ) . 

We now derive the channel — e. g. completely positive trace-preserving map — T which optimally approximates the 
(anti-linear) multi-phase conjugation on equatorial states, namely 



I' '<{"/!» = ^(10) + e**|l) + e*"|2) + • • • ) — W ({&})) = ^(|0) ~ ■ " '' |'l ) + < '-|2H i. (.(>) 



The map T is covariant under the multi-phase U(l) x( - d ^ group, namely 

T(CT({^}) p C7 f ({^})) = C/*({^}) T(p) t/ T ({<^}). (7) 

In the above equation O* (O t ) denotes the complex conjugation (transposition) of the operator O with respect to 
the orthonormal basis @ kept as real. Among all completely positive trace-preserving (CPT) maps satisfying the 
covariance condition J7J), we single out those maps which maximize the fidelity between the output state and the 
ideally transformed state in Eq. (JHJ) 



Ft 



(|^*({^})),^(I^({^-})>(^C{^})I)) -Tr [|^C{^}))(V^({0i})l ^(IVXi^DXV'd^})!)]- (8) 



Following Ref. J6(, we solve the optimization problem under the covariance condition JJJ, using the positive operator 
R on H ® H defined as 

fl=(T<8>I)|l»«l|, (9) 

where X is the identity channel, and = 53i=o ^ s * ne rnaximally entangled vector on 7i ® Ti relative to the 
orthonormal basis B in Eq. The correspondence T <-> i? is one-to-one and can be inverted as 

T(p) = Tr H2 [l(g>p T R}. (10) 

In terms of the operator R the trace-preservation condition for T reads 

Tr Hl [R] = t (11) 
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and the covariance property of the map J7J becomes the invariance for R 

[J2,tT({M) 82 ]=0. (12) 
This, in turn, via the Schur Lemma implies the following block form for R 

R= R a - (13) 

a: equiv. classes 

The index a runs over the equivalence classes of irreducible one-dimensional representations of E^({0.,})® 2 (without 
loss of generality we suppressed the complex conjugation). The equivalence classes with the respective characters are 
listed in Table [Q 



Equivalence Classes 


Characters 


|00) 


1 


|H) 








|01),|10) 









TABLE I: Equivalence classes and respective characters of irreducible one-dimensional representations of U({(j>j}) 
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Noticing that IV'oKV'ol = |V , o)(V , o|*i an( l using Eqs. ©, 11(J|) and (|12|l . we can rewrite the fidelity in Eq. (JSJ) as 

^ = Tr[|^ }(^or R\. (14) 

From Eq. Q, one has |^o) <^o|® 2 = d~ 2 YUj,k.l 1*^1® I and the fidelity ® is Cqual t0 F = d ~ 2 
where rj"] fc = (j7|i? a |«fc). As argued in Refs. 0,0, the maximum F is then obtained when the off-diagonal terms of 



the operator R are positive and as large as possible, namely when the blocks R a are rank-one or in equations 

1 



R 



Af 



J2 Cl \ii)(ii\ + Cij(aiij\ij) + Pij\ji)){a* i:j {ij\ + 



i>3 



(15) 



where M = Ci + Ei>j c ij)/d ^ s a normalization constant, |ay | + 1/9^1 2 = 1, and Cj ^ 0, Cy ^ 0, V i, j. The fidelity 
then takes the form 



F = Tr[|Vo>(V>or fl] 



d 



d 2 E^ + E 



'/ E, <■■ ■ E, ,, <\, 



^ c, + ]T cy (Kf + |/%| 2 + 2KK-/3*.)) 



(16) 



1 Ei>j ■Cjj^iaijPtj, 



Now, the maximum of 3t(oiij0*j) is achieved when onj = /3jj, implying 
Qy = 2 -1 / 2 . The fidelity becomes 



F 



Ei>j C «J 



(e* c » + e 



i>j c ij 



1/2 for all Therefore we choose 



(17) 
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In order to maximize F we put c, = for all i. The optimal fidelity takes the following simple form 

F-\ (18) 
Then, we impose the trace-preservation condition, obtaining 

Tn[i?] = — E^M + b'>0*l) 

= J2biM(i\ + \j){j\) (19) 

= a, 

where 

hi = d-^ • (20) 

Since the projector |0)(0|, for example, appears in the sum l|19[) multiplied by ^a=a kiQ> the coefficients are constrained 
as follows 

d-l 

$>o = l, (21) 

i=l 

and, similarly, for the element one has 

3-1 d-l 

E & ^ + E ^ = L ( 22 ) 

i=0 i=i+l 

Rearranging the positive coefficients {&ij}i>j into a square matrix array, they define the lower-trianguar section of a 
square matrix. Such a matrix can be uniquely completed to a null-diagonal symmetric bistochastic matrix, that is, a 
symmetric null-diagonal matrix with non-negative entries, such that all its rows' and columns' entries sum up to one, 
namely all its rows and columns are probability distributions. 
Up to now the operator R is simplified as follows 

R=J2 b ^ + ^))m+m), (23) 

where the coefficients hj's, uniquely defining a map T achieving the optimal fidelity 2/d, are the entries of a null- 
diagonal symmetric bistochastic matrix. From Eq. (|2.'jf) and the reconstruction formula (|10|) . an optimal phase 
covariant transposition map can then be easily expressed in the Kraus form as follows 

t(p) =EMI^1 + b'X*!) p CN>0'l + \M\)- (24) 

i>j 

Notice that the constraint Il22[l over {hj} is indeed very strong: for qubits and qutrits it suffices to completely and 
univocally determine the map. In the case of qubits the only null-diagonal symmetric bistochastic matrix is 

{hi} = (jj , (25) 

and the optimal transposition map is the unitary transformation 

T 2 (p) = a x pa x = p*, (26) 

which clearly achieves F=l. For d = 3 there is again a unique choice for a null-diagonal symmetric bistochastic 
matrix, which is given by 

/ 1/2 1/2^ 

{hi}= 1/2 1/2 | . (27) 
V 1/2 1/2 
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For d 4 there exist many optimal maps. For example, for d — 4 we can parametrise the family of maps by varying 
two positive parameters p\ and P2 such that < Pi + P2 = P12 < 1 ; 

/ 

Pi 



{bij} = 



Pi P2 

1 - Pl2 
P2 1 - Pl2 
\l-Pl2 P2 Pi 



-Pl2\ 
P2 
Pi 





(28) 



III. UNIVERSAL TRANSPOSITION AND MULTI-PHASE CONJUGATION 

In Refs. 0, H it was shown that the fidelity of the optimal universal transposition map is equal to the fidelity of 
the optimal universal pure state estimation |ld lllj , that for a single input copy is given by 

F = dh> ^ 

which is always lower than the fidelity of optimal phase covariant transposition l|18|l . as expected. The equivalence 
between transposition and state estimation means that an optimal universal transposition can be achieved by optimally 
estimating the input state and then preparing the transposed state. In this sense the optimal universal transposition 
is a "classical" map. 

In contrast to the universal case, the phase covariant transposition map cannot be achieved by phase estima- 
tion/preparation. Indeed, the fidelity of optimal multi-phase estimation for a single input copy is given by 



which is always smaller that the optimal fidelity of the phase covariant transposition map (|18|1 . Hence, the optimal 
phase covariant transposition is a genuinely quantum channel. The situation is particularly striking in the case of 
qubits, where it is possible to perfectly transpose all equatorial states, while the phase can never be measured exactly 
with finite resources |13| . 



IV. CONVEX STRUCTURE AND PHYSICAL REALIZATION OF OPTIMAL MAPS 



In Section [H] we have shown that optimal multi-phase conjugation maps are in one-to-one correspondence with 
null-diagonal symmetric bistochastic (NSB) matrices, which form a convex set. On the other hand, every bistochastic 
matrix is a convex combination of permutation matrices — this is the content of the Birkhoff theorem jl4j . The null- 
diagonal and symmetry constraints, however, force the convex set of NSB matrices to be strictly contained into the 
convex polyhedron of bistochastic matrices. This fact causes the extremal NSB matrices to eventually lie strictly 
inside the set of bistochastic matrices, generally preventing them from being permutations. 

The geometrical study of the set of NSB matrices and its extremal points can shed some light on the unusual feature 
that there exist different "equally optimal" maps. The problem arises for dimension at least d = 4. In this case the 
decomposition of the matrix {b^} in Eq. I|28[) into extremal components is 



{bij} =pi 



(° 


f 














f 






(0 








1\ 


1 











+ P2 











1 


+ P3 








f 














1 


f 














f 













1 


0) 






1 





0) 




V 








0/ 



Pi J 



p(D + p 2 p(2) +p ,p(3) 



P3J 



(31) 



where pi,p2,p% ^ and pi +P2 +P3 = 1 • A natural question is now which optimal maps can be achieved with minimal 
resources. 

In order to physically realize a given CPT map £ , one needs to design a specific unitary interaction U and prepare 
an ancilla in a specific state, say |0)(0| o , in such a way that 

f(p)=Tr o [C7(p®|0)(0| o )C/ t ]. (32) 

This is always possible 0, 0] . The existence of equivalently optimal maps allows us to choose between realizations 
with either a smaller ancilla dimension, or a more flexible ancilla state-preparation. 
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More explicitly, for d — 4, we define three unitaries Ui, U% and E/3 on C 

J 32 J 10/ \ J 31 J 20/ 

where = + Each of them realizes an extremal optimal multi-phase conjugation map (corresponding to 

Pk = 1 in Eq. (|3I|I for a given k), namely 




Tt\p) =Y, p <i )T »f ir V = Tr ^ (P® l°X°U) ^l. ( 34 ) 



where |0)(0| a is a fixed qubit ancilla state. Notice that the ancilla must not necessarily be in a pure state, and the 
optimal map is equivalently achieved for diagonal mixed ancilla state a|0)(0| Q + /3|l)(l| a . By adding a control qutrit, 
we can now choose among any of the optimal maps using the controlled-unitary operator on C 4 <g) C 2 <S> C 3 

U = U 1 <E> |0)(0| + U 2 <8> |1)(1| + U 3 ® |2)(2|. (35) 

Any optimal multi-phase conjugation map can now be written as 

%{p) = Tr Q , b [U (p® |0)(0|„ <8> (T h ) C/ f ] (36) 

where <7(, is a generic density matrix on C 3 . By superimposing or mixing the three orthogonal states {|0), |1), |2)} of 
the qutrit we control the weights Pi,P2,P3 hi Eq. I|31l) via the d iago nal entries of ov 
Eqs. (|33|I - H36|1 can be generalized for higher even dimensions [lij, with 

& _i 
2 1 

C^fc = ^2 T k®2i®2j,2i®2j ® fe = 1, . . . , d — 1, 

t,j=0 
d— 1 

C/ = ^f/ fc ® (37) 
fc=i 

ri fe) (p) = Tr a [f/ fc (p®|0)(0| Q )C/t], 
T d (p) = Tr a , b [(7 (p® |0)(0|„ O a b ) [7 f ] 



where Uk's are unitary operators acting on C d <8> C d / 2 , U is a control-unitary operator on <C d ® C d / 2 ® C d_1 , |0)(0| a is 
a fixed (d/2)-dimcnsional pure state, and ct;, is a generic (d— l)-dimensional density matrix. The minimum dimension 
of the ancilla space required to unitarily realize an optimal phase covariant transposition map is d/2, generalizing 
the result for d = 4, for which just a qubit is needed (see Eq. (|34H % 1. Notice that realization of phase covariant 
transposition generally needs much less resources than realization of universal transposition: the dimension d/2 of the 
ancilla space in the phase covariant case has to be compared with the dimension d 2 required in the universal case |lj . 

As a final remark, notice that unitary realizations of CPT channels are far from being uniquely determined |15j : 
here we chose the controlled unitary structure because of its clear geometrical interpretation in connection with the 
convex structure of the polyhedron of optimal maps. 



V. CONCLUSIONS 



In this paper we have derived the channel which optimally approaches in fidelity the time-reversal of multi-phase 
equatorial states in arbitrary (finite) dimension. We show that, in contrast to the customary case of the Universal- 
NOT on qubits (or the universal conjugation in arbitrary dimension), the optimal phase covariant time-reversal for 
equatorial states is a nonclassical channel, which cannot be achieved via a measurement /preparation procedure. We 
have given unitary realizations of the optimal time-reversal channel with minimal ancillary dimension, exploiting the 
simplex structure of the optimal maps. The optimal channels are related to null-diagonal symmetric bistochastic 
matrices. For d 4 this gives a simplex structure of equivalently optimal maps. 
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